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We present a theoretical analysis of excitation modes in Bose-Einstein condensates in ultracold 
alkali-metal gases for large scattering lengths and momenta where corrections to the mean field 
approximation become important. We assume that the effective interaction in the metastable, single 
channel, gaseous phase has a well defined Fourier transform that scales with the scattering length. 
Based on this we show that for increasing scattering lengths or equivalently increasing densities the 
system becomes less correlated and that at large values of the scattering length Bragg scattering 
measures directly the Fourier transform of the effective two-body potential. We construct model 
potentials which fit the recently measured line shifts in S5 Rb by Papp et al. (Phys. Rev. Lett. 
101, 135301 (2008)), and show that they fix the low momentum expansion of the effective range 
function. We find excellent agreement with the experimental data when the effective range is ^ 1 
and the coefficient of the fc 4 -term is — 7.5 ± 0.5 in scattering length units. The resolution in Bragg 
scattering experiments so far does not reveal details of the frequency dependence in the dynamic 
structure function S(k,u>) and we show that the Feynman spectrum determines the measured line 
shifts. We propose the possibility of a transition to a novel density wave state. 
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Bragg spectroscopy measurements of excitation spec- 
tra in ultracold atomic gases with Bose-Einstein con- 
densation have shown excellent agreement with the pre- 
dictions of the Bogoliubov theory of elementary excita- 
tionsi^i&i However, in recent experiments 5 - on 85 Rb 
particle density, scattering length and transferred mo- 
mentum have reached such high values that deviations 
from that are clearly observed and a possible relation to 
the 4 He-like roton behavior has been speculated. 

From the theoretical side the dynamic structure func- 
tion S(k, lj) determines the excitation modes of a quan- 
tum Bose system. At low momenta and frequencies typ- 
ical Bose systems show a sharp peak which defines the 
elementary excitation mode. At given k but higher fre- 
quency, multiphonon excitations form a broad continuous 
distribution. In the line shift experiments 5 -, the resolu- 
tion is not high enough to separate out these contribu- 
tions. Consequently, the experiments reveal information 
only on the average structure. For a fixed k, S(k,u>) sat- 
isfies strict sum rules 
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as the deviation of the maximum of the broad frequency 
distribution of the measured S(k,uj) from the free parti- 
cle value h 2 k 2 /(2m). If we identify the maximum with 
the average value (w(k)}, the line shift can be calculated 
from purely ground state quantities 
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with S(k) the static structure factor. These sum rules 
implies that the average value of the frequency (u) = 
mi(q) I (hma(q)) — hk 2 /(2mS(k)) obeys the Feynman 
spectrum. In the analysis of experiments 5 the frequency 
dependence of S(k,w) was fitted with a Gaussian func- 
tion. There, line shift at a given momentum k is obtained 



In this Letter we present results on the static and dy- 
namic structure functions of bosonic fluids using the cor- 
related basis function method 6 widely applied to strongly 
correlated fluids like the charged Bose gasi and in par- 
ticular superfluid 4 He 8 £. We assume that the effective 
interaction scales with the scattering length a by set- 
ting the length unit equal to a and the energy unit equal 
to h 2 /(2ma 2 ). Since we compare our results with mea- 
surements on 85 Rb, we set m equal to the 85 Rb mass. 
Experiments are done for the fixed momentum trans- 
fer k — 47r/780nm and density p — 7.6 x 10 13 cm~ 3 , 
but changing the scattering length a on a wide range 
153as < a < 891as where as is the Bohr radius. The 
only quantity then needed to control the many-body 
problem is the gas parameter x = pa 3 . 

Two simple models, which scale with the scattering 
length, have been studied extensively in the literature. 
The first one is the Fermi pseudopotential that acting on 
non-singular wave functions has the form (A-Kah 2 /m)S(r) 
with the strength proportional to the scattering length. 
The second one is the gas of hard spheres where the di- 
ameter of the particles is equal to the scattering length. 
These models give the linear behavior for the line shift 
as a function of the scattering length when k<l, but 
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FIG. 1: Comparison of the CBF S(k,u>) (black line), the CBF 
S(k, uj) folded with the instrumental resolution functions of 
Ref . [H (red line) , and the Feynman approximation also folded 
with these functions (blue line) for soft spheres potential with 
the radius R = 3.5a. Notice that the curves have been artifi- 
cially shifted for clarity. 



begin to deviate from each other and also from the ex- 
periments when ka > 0.4. That is why we extend in this 
work the analysis to more realistic effective potentials, 
which have finite range and finite strength, but assume 
that they have a Fourier transform. 

The analysis of the line shift experiments requires a 
theoretical model for S(k, u) that allows finding the spec- 
trum of elementary excitations. A suitable choice for 
that is provided by the Correlated Basis Function (CBF) 
method 9 , which is very well suited to describe weakly in- 
teracting Bose gasesi and yields reasonable results also 
for 4 He. We have found, however, that due to poor in- 
strumental resolution, the detailed structure of S(k, oS) 
is washed out in the experiments and that the elemen- 
tary excitation spectrum is very well described by the 
Feynman approximation. This is seen in Fig. [T] which 
shows a CBF calculation for a system of bosons inter- 
acting through the soft spheres potential V(r < R) = 
Vo, V(r > R) = for R = 3.5 in scattering length units. 
We have found that this simple model describes fairly 
well the experimental line shifts 5 . The figure also shows 
the CBF S{k 1 Lo) and the Feynman approximation, both 
folded with the instrumental resolution function of the 
experiments. The maxima of these broadened distribu- 
tions agree very well and we conclude that the relevant 
quantity describing the dynamics in the experiment is 
S{k). 

We model the gas of N 85 Rb atoms of mass m with 
the Hamiltonian 
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and the ground state many-body wave function 

N 

* (r u r 2 ,...,r N ) = J[ /fl^-i^l) (4) 

i=i<j 

is defined as a product of two-body correlation 
functions^, The unknown correlation function f(r) is the 
solution of the variational problem which minimizes the 
expectation value of the Hamiltonian. With diagram- 
matic techniques the resulting Euler-Lagrange equation 
can be written in a Bogoliubov-like form for the structure 
function S{k)~, which in scattering length units becomes 
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The effective, many-body potential in coordinate space 
splits into two parts, 



V p - h (r) = V(r) + 
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+ [g(r)-l]K d (r) + y(r)] + 2(vVsW) . 

The potential on the first line is just the bare two-body 
potential. In momentum space 
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defines the uniform limit approximation (with V p -h{k) = 
V(k))^^- for S(k), that reduces in turn to the ordi- 
nary Bogoliubov expression when V(r) is further approx- 
imated by the pseudopotential. The rest of the terms in 
Eq. ^ contain the many-body contributions to the ef- 
fective interaction. The radial distribution function g(r) 
is calculated from the inverse Fourier transform 
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and the induced potential w- ln d{r) is obtained by taking 
the inverse Fourier transform of the expression 



w ind (k) = -±(2S(k) + l) U-^y 



(9) 



Eqs. (UJ)-© form a set of self-consistent equations that 
have to be solved iteratively. 



A measure of the strength of the correlations in a quan- 
tum system is given by the amount g(r) deviates from 
unity. From Eq. ^ we see that many-body effects di- 
minish when g(r) — > 1. That happens when x increases 
for any potential with a Fourier transform as shown by 
Eqs. ijHJ) and provided that S(k) has no singulari- 
ties. For the PTGT-potential we find V p - h (k) « V(k) 
for x > 0.02. 
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FIG. 2: The line shift Aui(k) as a function of the scatter- 
ing length calculated using the PTGT (solid line) and SW 
(dashed line) potentials. Theoretical results are compared 
with the experimental points with error bars of Ref. H Also 
shown are the hard sphere (dashed-dotted curve) and the 
pseudopotential (dotted curve) results. The upper axis shows 
ka for a momentum 47r/780nm -1 . 



In the long wave length limit 



V p -h(k = 0) = mc 2 



(10) 



where c is the speed of sound. At small x, i.e. x < 10~ 3 , 
one recovers the pseudopotential result V^_h(0) = 87rx. 
The agreement with Monte Carlo simulations^ is very 
good up to x k 0.1 even for the hard spheres^ 2 -. For soft- 
core type potentials, the agreement should be even better 
at large x, because many-body contributions to V p -h(k) 
become less important. 

We use two models for the single channel potential. 
The first one is the Poschl- Teller potential to which we 
have added a weak, attractive, Gaussian tail (PTGT). 
The attractive tail mimics the long range part of the 
van der Waals interaction, and it improves the fit to the 
experimental line shifts, 
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where G x = 0.306038, G 2 = 0.012955, r = 7.377861, 
(i = 0.383580, and a = 4.799// 2 . The second model in 
our analysis is the step-well potential (SW) 
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with Vi = 0.110786, V 2 = 0.010572, R x = 4.684541 and 
i?2 = 1.5i?i. All parameters are in scattering length 
units. They are fixed by setting the effective range 
r e ff/a <C 1, based on the multichannel scattering anal- 
ysis near the Feshbach resonance by Bruun et alJ^ and 
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FIG. 3: In the upper panel the effective many-body poten- 
tials as a function of momentum are shown for the PTGT- 
potential and different values of the scattering ength ranging 
from 100a_g to 1500a s with 100as steps. The lowest curve 
is for the lowest value of a and the highest curve is equal to 
the Fourier transform of the bare two-body potential V(k). 
The lower panel shows the radial distribution functions cal- 
culated using the PTGT-potential for the same values of the 
scattering length a. The lowest curve is for the lowest value 
of a. 



by fitting the line shift data. 

In Fig. [2]we compare our results to experiments. Up to 
a = 400ao the behavior is almost linear and also both the 
pseudopotential and the hard sphere models fit well with 
the experimental data. Those models, though, are unable 
to reproduce the strong downwards bending observed in 
the line shifts at higher values of the scattering length. 

The effective interaction V p ~h(k) determines the line 
shifts completely within the resolution of the experi- 
ments. In Fig. we show its behavior for different 
values of the scattering length ranging from lOOas to 
1500a_B with 100a_B steps using the PTGT-potential of 
Eq. (fTTj) . Since we have fixed the density to the experi- 
mental value, that range corresponds to 10~ 5 < x < 0.04 
in the gas parameter. In the figure we have divided out 
the Bogoliubov constant 8irx and that is why at k = the 
lowest curve is very close to unity. Increasing the scatter- 
ing length increases the strength of the potential at k = 
since the speed of sound increases (see Eq. (ITU1) ). The 
upper curves at the largest values of a converge nicely to 
the Fourier transform of the bare potential. This means 
that many-body effects in V p -h{k) vanish and therefore 
the largest values available in experiments directly test 
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FIG. 4: The effective range function as a function of the two- 
body scattering momentum for the PTGT (diamonds) and 
SW (circles) potentials. Least square fits are done for the 
range < k < 0.5 to the PTGT (solid curve) and SW (dashed 
curve) potentials. Also shown are the results for the hard 
sphere (dashed-dotted curve) and <5-function (dotted curve) 
potentials. 

the Fourier transform of the effective, bare potential. As 
long as the system remains in the gaseous phase, the 
particle-hole potential is directly given by the bare po- 
tential even if one keeps increasing a. 

The momentum transfer fixed in the experiments 
varies in 1/a-units from 0.13 to 0.76 as depicted in Fig. 
[3J Clearly understanding the behavior of the effective 
interaction as a function of the momentum transfer is es- 
sential since at ka ~ 1 the effective potential approaches 
zero defining an inverse healing distance £ _1 where the 
line shift should also go to zero. 

The potential models chosen in Eqs. (fTTj) and (fT2")l give 
very similar line shifts within the experimental range, but 
at larger a they behave very differently. Oscillations ev- 
ident in the Fourier transform of the SW potential make 
it go through zero at ka w 0.95 whereas the PTGT- 
potential remains positive up to ka < 1.7. That dif- 
ference could be experimentally tested by increasing the 
used momentum transfer by ~20%. The SW potential 
also leads to the density wave instability at a w 2260as 
with ka = 1.15 whereas with the PTGT-potential the 
system remains in the gaseous phase. It could also be 
interesting to search for that phase transition experimen- 
tally. 

In order to better understand the importance of the 
correlations we have plotted in Fig. [3J3 the radial distri- 
bution function g(r) of the PTGT-potential for scatter- 
ing lengths spanning the range lOOas < a < 1500as. It 
shows that correlations diminish with increasing scatter- 
ing length. That is why the effective potential Vp_h(fc) is 
completely determined by the bare two-body potential at 
large a. In some sense that is equivalent to the behavior 
of the charged Bose gas^ where the high density limit is 
the weakly correlated fluid. 



In comparison of different models we plot in Fig. [4] the 
s-wave effective range function k cot((5o (k j) defined by the 
on-shell T-matrix for the different models used. At small 
k < 0.25 it is very close to -1 supporting the scattering 
length approximation. At larger k our model potentials 
give very similar downward bending. At k > 0.6 the first 
order Born approximation becomes very accurate and we 
can conclude that the Fourier transform of the bare two- 
body potential V(k), determines completely the large k 
behavior of the effective range function. On the other 
hand, V p -h(k) for these momenta is also accurately de- 
termined by V(k). This is not surprising since at large 
momenta the physics is governed by the short-distance 
two-body problem. This means that at large scattering 
lengths V(k) could be determined by the line shift mea- 
surements. In Fig. Q]we also plot the hard sphere poten- 
tial result and as shown it bends to the wrong direction. 
One can understand the link between V p -h(k) and the 
T-matrix by realizing that the large k behaviour of both 
quantities is determined by the 2-body physics at short 
distances 

Finally we give a simple parametrization by making a 
least square fit to the calculated effective range functions 
with a fourth degree polynomial 

kcot{S (k)) = -l + ^r cS k 2 + Pk 4 (13) 

when < k < 0.5. The k A coefficient is the shape param- 
eter P, which turns out to be important for reproducing 
well the bending of the line shift at large a. With the 
PTGT-potential we find P = -7.28 and with the SW- 
potential P = —7.88. The quality of the polynomial fits 
is also show in Fig. [5] 

In summary, we have assumed that the effective two- 
body potential scales with the scattering length and has a 
Fourier transform. Within these assumptions many-body 
correlations diminish with increasing scattering length 
and line shift experiments could be used to measure di- 
rectly the Fourier transform of the bare two-body ef- 
fective potential. The downward bending of the line 
shift clearly shows that effective interactions with finite 
strength and range are needed. By slightly increasing the 
gas parameter and momentum transfer one could study 
experimentally the formation of a possible density wave 
instability and by increasing the experimental resolution 
one could separate the elementary excitation mode from 
the continuum of multiphonon contributions as seen in 
4 He and the charged Bose gas. 
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